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Abstract 

We consider the Ginzburg-Landau MTV-model that describes M N- 
vector cubic models with 0(M)-symmetric couplings. We compute the 
renormalization-group functions to six- loop order in d = 3. We focus on 
the limit N — » which describes the critical behaviour of an M-vector model 
in the presence of weak quenched disorder. We perform a detailed analysis 
of the perturbative series for the random Ising model (M = 1). We ob- 
tain for the critical exponents: 7 = 1.330(17), v = 0.678(10), n = 0.030(3), 
a = -0.034(30), (3 = 0.349(5), u = 0.25(10). For M > 2 we show that 
the 0(M) fixed point is stable, in agreement with general non-perturbative 
arguments, and that no random fixed point exists. 

PACS Numbers: 75.10.Nr, 75.10.Hk, 02.30.Lt, 64.60.Ak, ll.10.-z, 64.60.Fr, 
75.40.Cx. 
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I. INTRODUCTION. 



The critical behavior of systems with quenched disorder is of considerable theoretical 
and experimental interest. A typical example is obtained by mixing an (anti)-ferromagnetic 
material with a non-magnetic one, obtaining the so-called dilute magnets. These materials 
are usually described in terms of a lattice short-range Hamiltonian of the form 

?4 = ~ J Pi Pi 8* ' ( L1 ) 

<ij> 

where Sj is an M-component spin and the sum is extended over all nearest-neighbor sites. 
The quantities pi are uncorrelated random variables, which are equal to one with probability 
x (the spin concentration) and zero with probability 1 — x (the impurity concentration). The 
pure system corresponds to x — 1. One considers quenched disorder, since the relaxation 
time associated to the diffusion of the impurities is much larger than all other typical time 
scales, so that, for all practical purposes, one can consider the position of the impurities 
fixed. For sufficiently low spin dilution 1 X. 1.6. clS long as one is above the percolation 
threshold of the magnetic atoms, the system described by the Hamiltonian 7i x undergoes a 
second-order phase transition at T c (x) < T c (x = 1) (see e.g. Ref. [I] for a review). 

The relevant question in the study of this class of systems is the effect of the disorder 
on the critical behavior. The Harris criterion states that the addition of impurities 
to a system which undergoes a second-order phase transition does not change the critical 
behavior if the specific-heat critical exponent a pure of the pure system is negative. If a pure 
is positive, the transition is altered. Indeed, in disordered systems the exponent v satisfies 
the inequality v > 2/d [|3|,|3| - this fact has been questioned however in Refs. |5]|| - and 
therefore, by hyperscaling, a ran dom is negative. Thus, if a pur c is positive, a ran dom differs 
from a pure , so that the pure system and the dilute one have a different critical behavior. In 
pure M- vector models with M > 1, the specific-heat exponent a purc is negative; therefore, 
according to the Harris criterion, no change in the critical asymptotic behavior is expected 
in the presence of weak quenched disorder. This means that in these systems disorder leads 
only to irrelevant scaling corrections. Three-dimensional Ising systems are more interesting, 
since a pure is positive. In this case, the presence of quenched impurities leads to a new 
universality class. 

Theoretical investigations, using approaches based on the renormalization group ]?|-|57|, 
and numerical Monte Carlo simulations |38|-|52"f , support the existence of a new random Ising 
fixed point describing the critical behavior along the T c (x) line: the critical exponents are 
dilution independent (for sufficiently low dilution) and different from those of the pure Ising 
model. 

Experiments confirm this picture. Cristalline mixtures of an Ising-like uniaxial antifer- 
romagnet with short-range interactions (e.g. FeF2, MnF 2 ) with a nonmagnetic material 
(e.g. ZnF 2 ) provide a typical realization of the random Ising model (RIM) (see e.g. Refs. 
|53|-|69|). Some experimental results are reported in Table |. This is not a complete list, 
but it gives an overview of the experimental state of the art. Other experimental results 
can be found in Refs. flU36fl. The experimental estimates are definitely different from the 



values of the critical exponents for pure Ising systems, which are (see Ref. |70| and references 
therein) 7 = 1.2371(4), v = 0.63002(23), a = 0.1099(7), and (3 = 0.32648(18). Moreover, 
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TABLE I. Experimental estimates of the critical exponents for systems in the RIM universality 



class. 


Ref. 


materia,! 


concentration 


/ 






3 




r^/ni^r^ 


x = U.4o 










@ 


Mn^Zn^Fs 


x = 0.75 


1.364(76) 


0.715(35) 






@ 


Fe 3; Zni„ :c F2 


x = 0.9 








0.350(9) 


m 


Mn^Zn^Fs 


x = 0.40,0.55,0.83 






-0.09(3) 






Mn a; Zni_ a; F2 


x = 0.5 








0.33(2) 




Fe x Zni_ :c F2 


x = 0.5 








0.36(2) 




Fe I Zni_ x F 2 


x = 0.93 






-0.10(2) 




n 


Fe a: Zni_ x E2 


x = 0.93 


1.34(6) 


0.70(2) 







they appear to be independent of the concentration. We mention that in the presence of an 
external magnetic field along the uniaxial direction, dilute Ising systems present a different 
critical behavior, equivalent to that of the random-field Ising model |7T[]. This is also the 
object of intensive theoretical and experimental investigations (see e.g. Refs. JF^,[F^]). 

Several experiments also tested the effect of disorder on the A-transition of 4 He that 
belongs to the XY universality class, corresponding to M = 2 |]7^-f79|]. They studied the 
critical behaviour of 4 He completely filling the pores of porous gold or Vycor glass. The 
results indicate that the transition is in the same universality class of the A-transition of the 
pure system in agreement with the Harris criterion [|80|. 



The starting point of the field-theoretic approach to the study of ferromagnets in the 



presence of quenched disorder is the Ginzburg-Landau- Wilson Hamiltonian [12 



H 



;i.2) 



where r cc T — T c , and ip(x) is a spatially uncorrelated random field with Gaussian distri- 
bution 



1 



exp 



4w 



;i.s) 



We consider quenched disorder. Therefore, in order to obtain the free energy of the system, 
we must compute the partition function Z(ip,g ) for a given distribution ip(x), and then 
average the corresponding free energy over all distributions with probability P(ip). By using 
the standard replica trick, it is possible to replace the quenched average with an annealed 
one. First, the system is replaced by N non-interacting copies with annealed disorder. Then, 
integrating over the disorder, one obtains the Hamiltonian |T2 



Ti 



AIN 



d d x 



+ r< 



9 ■ 

a,i 



ij,ab 




;i.4) 



where a, b = 1, ...M and i, j = 1, ...N. The original system, i.e. the dilute M-vector model, 
is recovered in the limit iV — >• 0. Note that the coupling w is negative (being proportional 
to minus the variance of the quenched disorder), while the coupling v is positive. 
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In this formulation, the critical properties of the dilute M-vector model can be inves- 
tigated by studying the renormalization-group flow of the Hamiltonian (|1.4j ) in the limit 
iV — > 0, i.e. of TCmo- One can then apply conventional computational schemes, such as 
the e-expansion, the fixed-dimension d = 3 expansion, the scaling-field method, etc... In 
the renormalization-group approach, if the fixed point corresponding to the pure model is 
unstable and the renormalization-group flow moves towards a new random fixed point, then 
the random system has a different critical behavior. 

It is important to note that in the renormalization-group approach one assumes that 
the replica symmetry is not broken. In recent years, however, this picture has been ques- 
tioned p6|-p8| on the ground that the renormalization-group approach does note take into 
account other local minimum configurations of the random Hamiltonian (|1.2p , which may 
cause the spontaneous breaking of the replica symmetry. In this paper we assume the valid- 
ity of the standard renormalization-group approach, and simply consider the Hamiltonian 
(fD§ for N = 0. 

For generic values of M and N, the Hamiltonian Hmn describes M coupled N- vector 



models and it is usually called MN model [TBI]. Hmn has four fixed points: the trivial 



Gaussian one, the 0(M)-symmetric fixed point describing N decoupled M-vector models, 
the 0(MiV)-symmetric and the mixed fixed point. The Gaussian one is never stable. The 



stability of the other fixed points depends on the values of M and N (see e.g. Ref. [13] for 
a discussion). The stability properties of the decoupled O(M) fixed point can be inferred 
by observing that the crossover exponent associated with the 0(MiV)-symmetric interac- 
tion (with coupling it ) is related to the specific- heat critical exponent of the O(M) fixed 
point |^,[nj. Indeed, at the 0(M)-symmetric fixed point one may interpret Hmn as the 
Hamiltonian of iV M-vector systems coupled by the 0(MiV)-symmetric term. But this 
interaction is the sum of the products of the energy operators of the different M-vector 
models. Therefore, at the O(M) fixed point, the crossover exponent associated to the 
0(MiV)-symmetric quartic term should be given by the specific- heat critical exponent (%m 
of the M-vector model, independently of N. This argument implies that for M = 1 (Ising- 
like systems) the pure Ising fixed point is unstable since <ft — a i > 0, while for M > 1 the 
O(M) fixed point is stable given that cxm < 0. This is a general result that should hold 
independently of N. 

For the quenched disordered systems described by the Hamiltonian Hmo, the physically 
relevant region for the renormalization-group flow corresponds to negative values of the 
coupling u (T^HIl • Therefore, for M > 1 the renormalization group flow is driven towards the 
pure O(M) fixed point, and the quenched disorder yields correction to scaling proportional 
to the spin dilution and to |t| Ar with A r = — cum- Note that for the physically interesting 
two- and three- vector models the absolute value of cum is ver Y small: a 2 — —0.013 (see e.g. 
the recent results of Refs. [PSSH]) and "3 ^ -0.12 (see e.g. SH)- Thus disorder 
gives rise to very slowly-decaying scaling corrections. For Ising-like systems, the pure Ising 
fixed point is instead unstable, and the flow for negative values of the quartic coupling u 
leads to the stable mixed or random fixed point which is located in the region of negative 
values of u. The above picture emerges clearly in the framework of the e-expansion, although 
for the Ising-like systems the RIM fixed point is of order y/e ]TO] rather than e. 

The other fixed points of the Hamiltonian Hmo are located in the unphysical region 
u > 0. Thus, they are not of interest for the critical behavior of randomly dilute spin 
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models. For the sake of completeness, we mention that for M > 1 the mixed fixed point is 
in the region of positive u and is unstable [[HJ. The last fixed point is on the positive v — 
axis, is stable and corresponds to the (MiV)-vector theory for N — > 0. It is therefore in the 
same universality class of the self-avoiding walk model (SAW). Figure [1| sketches the flow 
diagram for Ising (M = 1) and multicomponent (M > 1) systems. 

The Hamiltonian TLmn has been the object of several field-theoretic studies, especially 
for M = 1, the case that describes the RIM. Several computations have been done in the 
framework of the e-expansion |J3j and of the fixed-dimension d = 3 expansion [93]. In these 
approaches, since field-theoretic perturbative expansions are asymptotic, the resummation 
of the series is essential to obtain accurate estimates of physical quantities. For pure systems 
described by the Ginzburg-Landau- Wilson Hamiltonian one exploits the Borel summabil- 
ity [9J| of the fixed-dimension expansion (for which Borel summability is proved) and of the 
e-expansion (for which Borel summability is conjectured), and the knowledge of the large- 
order behavior of the series |96| , |97[| . Resummation procedures using these properties lead to 
accurate estimates (see e.g. Refs. p8|- p.00|j90|j9^] ). 

Much less is known for the quenched disordered models described by Hmo- Indeed, the 
analytic structure of the corresponding field theory is much more complicated. The zero- 
dimensional model has been investigated in Ref. ||101|| . They analyze the large-order behavior 
of the double expansion in the quartic couplings u and v of the free energy and show that 
the expansion in powers of v, keeping the ratio A = —u/v fixed, is not Borel summable. 



In Ref. [[L 02 ], it is shown that the non-Borel summability is a consequence of the fact that, 
because of the quenched average, there are additional singularities corresponding to the ze- 
roes of the partition function Z(ip,g Q ) obtained from the Hamiltonian fll.2|) . Recently the 
problem has been reconsidered in Ref. [ |103| |. In the same context of the zero-dimensional 
model, it has been shown that a more elaborate resummation gives the correct determina- 
tion of the free energy from its perturbative expansion. The procedure is still based on a 
Borel summation, which is performed in two steps: first, one resums in the coupling v each 
coefficient of the series in u; then, one resums the resulting series in the coupling u. There 
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is no proof that this procedure works also in higher dimensions, since the method relies on 
the fact that the zeroes of the partition function stay away from the real values of v. This 
is far from obvious in higher- dimensional systems. 

At present, the most precise field-theoretic results have been obtained using the fixed- 
dimension expansion in d — 3. Several quantities have been computed: the critical expo- 
nents 



18.20.23-25.27.32.34.35,37 



R{ 



the equation of state |28| and the hyper universal ratio 
3"3"[ . The most precise estimates of the critical exponents for the RIM have been ob- 
tained from the analysis of the five-loop fixed- dimension expansion, using Pade-Borel-Leroy 
approximants |J7| . In spite of the fact that the series considered are not Borel summable, 
the results for the critical exponents are stable: they do not depend on the order of the 
series, the details of the analysis, and, as we shall see, are in substantial agreement with our 
results obtained following the precedure proposed in Ref. ||103|| . This fact may be explained 
by the observation of Ref. [|101|| that the Borel resummation applied in the standard way 
(i.e. at fixed v/u) gives a reasonably accurate result for small disorder if one truncates the 
expansion at an appropriate point, i.e. for not too long series. 

The MN model has also been extensively studied in the context of the e-expansion 
[|T|Jg| p~6| , p~9| , |2T] , |22| , |26| , |30|] . The critical exponents have been computed to three loops for generic 

Several studies also considered the 
In spite of these 
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values of M, N |22[ and to five loops for M 

equation of state [n,19,26[] and the two-point correlation function fH],21 



efforts, studies based on the e-expansion have been not able to go beyond a qualitative de- 
scription of the physics of three-dimensional randomly dilute spin models. The ^-expansion 
T0|1 turns out not to be effective for a quantitative study of the RIM (see e.g. the anal- 



ysis of the five-loop series done in Ref. JJtJ). A strictly related scheme is the so-called 
minimal-subtraction renormalization scheme without e-expansion [|105|| . The three- loop [29 
and four-loop |^T].fJl|.fJB| results are in reasonable agreement with the estimates obtained 
by other methods. At five loops, however, no random fixed point can be found |]36| using 
this method. This negative result has been interpreted as a consequence of the non-Borel 
summability of the perturbative expansion. In this case, the four-loop series could represent 
the "optimal" truncation. We also mention that the Hamiltonian (1A) for M — 1 and 
N — > has been studied by the scaling-field method 

The randomly dilute Ising model (|1.1|) has been investigated by many numerical simula- 

The first simulations were apparently finding critical exponents 

79| that this could be simply 



tions (see e.g. Refs. p8|-[51] 



depending on the spin concentration. It was later remarked 
a crossover effect: the simulations were not probing the critical region and were comput- 
ing effective exponents strongly biased by the corrections to scaling. Recently, the critical 
exponents have been computed |)TJ using finite-size scaling techniques. They found very 
strong corrections to scaling, decaying with a rather small exponent uo ~ 0.37(6), — corre- 
spondingly A = lov = 0.25(4) — which is approximately a factor of two smaller than the 
corresponding pure-case exponent. By taking into proper account the confluent corrections, 
they were able to show that the critical exponents are universal with respect to variations of 
the spin concentration in a wide interval above the percolation point. Their final estimates 
are reported in Table 0. 

In this paper we compute the renormalizat ion-group functions of the generic MN model 
to six loops in the framework of the fixed- dimension d = 3 expansion. We extend the 
three-loop series of Ref. p3 and the expansions for the cubic model (M = 1) reported 
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TABLE II. Best theoretical estimates of the critical exponents for the RIM universality class. 
We report results for the massive scheme in fixed dimension d = 3 ("d = 3 exp."), the minimal 
subtraction scheme without e-expansion ("e? = 3 MS"), and the best Monte Carlo ("MC") results. 









Method 


7 


V 


V 


UJ 


This 1 
Ref. |7 
Ref. [34,3 
Ref. |l 


vork 
|, 2000 
6|, 1999 
|, 1998 


d = 
d = 
d = 


3 exp. 0(g 6 ) 
3 exp. 0{g 5 ) 
-- 3 MS 0(g 4 ) 
MC 


1.330(17) 
1.325(3) 

1.318 
1.342(10) 


0.678(10) 
0.671(5) 
0.675 
0.6837(53) 


0.030(3) 
0.025(10) 

0.049 
0.0374(45) 


0.25(10) 
0.32(6) 
0.39(4) 
0.37(6) 



in Ref. |37| (five loops) and Ref. |92| (six loops). We will focus here on the case N = 
corresponding to disordered dilute systems. Higher values of N are of interest for several 
types of magnetic and structural phase transitions and will be discussed in a separate paper. 



For M = 1 and N > 2 the six-loop series have already been analyzed in Ref. [92| where we 
investigated the stability of the 0(iV)-symmetric point in the presence of cubic interactions. 
We should mention that two-loop and three-loop series for the MN model in the fixed 
dimension expansion for generic values of d have been reported in Refs. p7| , |32| . 

For M — 1, N — 0, we have performed several analyses of the perturbative series follow- 
ing the method proposed in Ref. ||103[ |. The analysis of the /3-functions for the determinaton 
of the fixed point is particularly delicate and we have not been able to obtain a very robust 
estimate of the random fixed point. Nonetheless, we derive quite accurate estimates of the 
critical exponents. Indeed, their expansions are well behaved and largely insensitive to the 
exact position of the fixed point. Our final estimates are reported in Table [IT], together with 
estimates obtained by other approaches. The errors we quote are quite conservative and are 
related to the variation of the estimates with the different analyses performed. The over- 
all agreement is good: the perturbative method appears to have a good predictive power, 
in spite of the complicated analytic structure of the Borel transform that does not allow 
the direct application of the resummation methods used successfully in pure systems. For 
M > 2 and iV = we have verified that no fixed point exists in the region u < and that 
the 0(M)-symmetric fixed point is stable, confirming the general arguments given above. 

The paper is organized as follows. In Sec. [H| we derive the perturbative series for 
the renormalization-group functions at six loops and discuss the singularities of the Borel 
transform. The results of the analyses are presented in Sec. [TTT] and the final numerical 
values are reported in Table |Tj. 



II. THE FIXED-DIMENSION PERTURBATIVE EXPANSION OF THE 
THREE-DIMENSIONAL MN MODEL. 

A. Renormalization of the theory. 

The fixed-dimension 4 field-theoretic approach |94| provides an accurate description 



of the critical properties of 0(iV)-symmetric models in the high-temperature phase (see 
e.g. Ref. ||100|| ). The method can also be extended to two-parameter 4 models, such as 



the MN model. The idea is to perform an expansion in powers of appropriately defined 
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zero-momentum quartic couplings. The theory is renormalized by introducing a set of zero- 
momentum conditions for the (one-particle irreducible) two-point and four-point correlation 
functions: 



r% J ( P ) = S ai ^ 1 \m 2 + p 2 + 0(p 4 



where 5. 



ai,bj 



SnhS, 



ab"ij j 



(2.1) 



(2.2) 



and 



S a i,bj,ck,dl — - (8 a i t bj8ck,dl + ^ai,ck^bj,dl + &ai,dl$bj ,ck) 



a 



ai,bj,ck,dl 



8ij8ik8il 77 {&ab&cd + &ac&bd + ^ad^t 



be) 



(2.3) 
(2.4) 



Eqs. Q2.1D and ( [2.2|) relate the second-moment mass m, and the zero-momentum quartic 
couplings u and v to the corresponding Hamiltonian parameters r, uq and vq\ 



u q = muZuZ^ 2 , v = mvZ v Z 4> 2 . 

In addition we define the function Z t through the relation 

^ai,bj(fy — $ai,bjZ t \ 



(2.5) 



(2.6) 



where r^ 1 ' 2 ) is the (one-particle irreducible) two-point function with an insertion of |0 2 . 

From the pertubative expansion of the correlation functions T^ 2 \ and I^ 1,2 ) and the 
above relations, one derives the functions Z^u, v), Z u (u,v), Z v (u,v), Z t (u,v) as a double 
expansion in u and v. 

The fixed points of the theory are given by the common zeros of the /3-functions 



(3 u (u,v) = m 
(3 v (u,v) 



m 



du 
dm 

dv 
dm 



(2.7) 



«();«() 



calculated keeping w and Vq fixed. The stability properties of the fixed points are controlled 
by the matrix 



( d(3 u (u,v) d/3 u (u,v) 
du dv 

8f3 v {u,v) d(3 v {u,v) 
du dv 



(2.8) 



computed at the given fixed point: a fixed point is stable if both eigenvalues are positive. 
The eigenvalues iOi are related to the leading scaling corrections, which vanish as £ _£Ji ~ \t 
where Aj = vu^. 

One also introduces the functions 
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dlnZ^ dlnZ^ dlnZ^ 

V<t>(u,v) = — - = /3 M a y , (2.9) 

a In m aii at> 

91n^ 91n^ 01nZ t 

a mm ou ov 



Finally, the critical exponents are obtained from 



r) = r]4(u*,v*), (2.11) 
v=[2-7I4(u*,v*)+t H {u*,v*)]- 1 , (2.12) 
1 = u(2-r l ). (2.13) 



B. The perturbative series to six loops. 

We have computed the perturbative expansion of the correlation functions ( |2.1p , ( |2.2| ) 
and (|2.6|) to six loops. The diagrams contributing to the two-point and four-point functions 



to six- loop order are reported in Ref. ||106|| : they are approximately one thousand, and it is 
therefore necessary to handle them with a symbolic manipulation program. For this purpose, 
we wrote a package in Mathematica ||107 |. It generates the diagrams using the algorithm 



described in Ref. [108|, and computes the symmetry and group factors of each of them. 
We did not calculate the integrals associated to each diagram, but we used the numerical 
results compiled in Ref. [ |106| ] . Summing all contributions we determined the renormalization 
constants and all renormalization-group functions. 

We report our results in terms of the rescaled couplings 

1()TT 167T 

u = — R MN u, v = — R M v, (2.14) 

where Rk = 9/(8 + K), so that the /3-functions associated to u and v have the form 
fc(u, 0) = -u + u 2 + 0(u 3 ) and fy(0, v) = -v + v 2 + 0(v 3 ). 
The resulting series are 

A = -u + + 2 ? + y W - 4 H 9 ° + Jtf (2.15) 



M 27(8 + MN) 

400(2 + M) _ 2 _ 92(2 + M) 



-2- Ji l iTlu ; 2 i — u(. u )-i-i 

uv — T7\2 UV ^~ u L^, "" 



27(8 + MN) (8 + M) 27(8 + M) 



u v 



i+j>3 



_ -,- 2| 12 4(190 + 41M)_ 3 16(131 + 25M) __ 2 

P * = ~ V + V + sTMN™ - 27(8 + M) 2 V ~ 27(8 + MN)(8 + M) UV (2 " 16) 

27(8 + MN) 2 ^ 3 iJ 

= 8 ( 2 ± MiV ) u 2 + 16(2 + M ^ uv + 8(2 + M) v 2 + y e<*W (2 17) 

^ 27(8 + MiV) 2 27(8 + MiV)(8 + M) 27(8 + M) 2 <4>3 
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2 + MN_ 2 + M_ 2(2 + MN) _ 2 

nn, — 11 — y -U 11 

' 8 + MN 8 + M (8 + MiV) 2 

4(2 + M) , 2(2 + M) a ^ „, 



'(8 + MiV)(8 + Mr + iw V2 + 53 e ^ V - (2 ' 18) 

For 3 < % + j < 6, The coefficients , 6^, and e-^ are reported in the Tables PTT| , [TV| , 
[V] and respectively. 

We have performed the following checks on our calculations: 

(i) Pu(u, 0), ^(m, 0) and r) t (u, 0) reproduce the corresponding functions of the O(MN)- 
symmetric model p8|J109|| ; 

(ii) /3^(0, v), 7)^(0, v) and r; t (0, v) reproduce the corresponding functions of the O(M)- 
symmetric model |p8|J109|| ; 

(iii) For M — 1, the functions /3„, 77^ and rjt reproduce the corresponding functions of 
the A^-component cubic model f9~2]| ; 

(iv) The following relations hold for N — 1: 

Az(it, x — w) + a; — it) = /^(0, 2), (2.19) 
^(m, x - u) = ?fa(0,:r), 
97t(rx, a; -u) = r) t (0,x). 



C. Borel summability and resummation of the series. 

Since field-theoretic perturbative expansions are asymptotic, the resummation of the 
series is essential to obtain accurate estimates of the physical quantities. 

In the case of the 0(iV)-symmetric 4 theory the expansion is performed in powers of the 
zero-momentum four-point coupling g. The large-order behavior of the series S(g) = s k9 k 
of any quantity is related to the singularity gi> of the Borel transform closest to the origin. 
Indeed, for large k, the coefficient behaves as 

s k ~ k\ (-a) k k h \l + 0(k~ l )\ with a = -l/g b . (2.20) 

The value of depends only on the Hamiltonian, while the exponent b depends on which 
Green's function is considered. The value of can be obtained from a steepest-descent 
calculation in which the relevant saddle point is a finite-energy solution (instanton) of the 
classical field equations with negative coupling p6[P7[ . Since the Borel transform is singular 
for g = gb, its expansion in powers of g converges only for \g\ < \gt\- An analytic extension 
can be obtained by a conformal mapping |9"9]] , such as 



> — g/db— 1 

y(g) = h — • (2.21) 
-g/9b + i 
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TABLE III. The coefficients b\f , cf. Eq.pl5|) 



i i B~ l R~ J b 



0.27385517 + 0.075364029 MN + 0.0018504016 Af 2 N 2 



0.4516155 + 0.22580775 Af + 0.018235606 MN + 0.0091178029 M 2 N 



0.25029007 + 0.15182596 M + 0.013340463 Af 2 + 0.0017061432 M N + 0.00085307161 Af 2 N 



0.051324521 + 0.03463704 Af + 0.00448739 Af 2 



-0.27925724 - 0.091833749 M N - 0.0054595646 Af 2 N 2 + 0.000023722893 Af 3 N' A 
-0.62922442 - 0.31461221 Af - 0.055501872 MN - 0.027750936 IP N 
+0.00041240116 Af 2 N 2 + 0.00020620058 Af 3 N 2 



-0.55267483 - 0.36698183 Af - 0.045322209 Af 2 - 0.0092541768 Af N 
-0.0036796747 M 2 N + 0.00047370687 Af 3 N 



-0.23360529 - 0.16634719 Af - 0.02383493 Af 2 + 0.00046867235 Af 3 
-0.0010494477 MN- 0.00065590478 M 2 N - 0.000065590478 Af 3 N 
-0.040934998 - 0.029976629 M - 0.0046082314 M 2 + 0.000073166787 M :i 



0.35174477 + 0.13242502 MN + 0.011322026 M 2 N 2 + 0.000054833719 Af a N* 

+8.6768933 10~ 7 M 4 N 4 

1.0139338 + 0.50696692 M + 0.12967024 M N + 0.064835121 M 2 N + 0.00073857427 M 2 N 2 

+0.00036928714 M 3 N 2 + 0.000021186521 Af 3 N s + 0.000010593261 M 4 N s 

1.2323698 + 0.85537154 M + 0.11959332 M 2 + 0.041310322 MN + 0.02290239 M 2 N 
+0.0011236145 M 3 N + 0.00017140951 Af 2 jV 2 + 0.00015195691 Af 3 N 2 + 0.000033126079 M 4 N 2 
0.80970369 + 0.61305431 M + 0.10724129 M 2 + 0.0015700281 Af 3 + 0.0058598808 MN 



+0.0041006979 M 2 N + 0.00067267459 Af 3 N + 0.000043647916 Af 4 N 

0.28991291 + 0.22599317 M + 0.04307445 Af 2 + 0.0013403006 Af 3 + 0.00003112665 Af 4 
+0.0007752054 M N + 0.0005229621 Af 2 N + 0.000055080316 Af 3 N - 6.2996919 10~ 6 Af 4 N 
0.044655379 + 0.035455913 M + 0.0071030237 Af 2 + 0.00027580942 Af 3 + 3.1767371 10~ b Af 4 
-0.51049889 - 0.21485252 MN - 0.023839375 M 2 N 2 - 0.00050021682 Af 3 N 6 

+2.0167763 10~ 6 Af 4 jV 4 + 4.4076733 10~ 8 Af 5 jV 5 ^ 

-1.7989389 - 0.89946945 Af - 0.30045501 Af Af - 0.15022751 Af 2 Af - 0.007214312 Af 2 N 2 
-0.003607156 M 3 Af 2 + 0.000038644454 M 3 Af 3 + 0.000019322227 Af 4 Af 3 + 1.3676971 10~ 6 Af 4 Af 4 

+6.8384853 10~ 7 Af 5 Af 4 

-2.8025568 - 2.0077582 Af - 0.3032399 Af 2 - 0.15673528 Af Af - 0.10009354 Af 2 Af 
-0.010862953 M 3 N - 0.00057173496 Af 2 Af 2 - 0.00013491592 Af 3 Af 2 + 0.000075475779 Af 4 Af 2 

+7.5116588 10~ 6 Af 3 Af 3 + 9.1864196 10~ 6 Af 4 Af 3 + 2.7152951 10~ 6 Af 5 Af 3 

-2.5110555 - 2.0026222 Af - 0.40321804 Af 2 - 0.014835404 Af 3 - 0.038290982 Af Af 
-0.026220763 M 2 N - 0.0032211897 Af 3 Af + 0.0001582232 Af 4 Af - 0.000012833973 Af 2 Af 2 

+0.000015558757 Af 3 Af 2 + 0.000020698159 Af 4 Af 2 + 4.8551437 10~ 6 Af 5 Af 2 

-1.3671792 - 1.1287738 Af - 0.24696793 M 2 - 0.011904064 M A + 0.00014192829 Af 4 
-0.0053723551 MAf - 0.0037885035 Af 2 Af - 0.00047438677 Af 3 Af + 0.000047179342 Af 4 Af 

+4.3956195 !Q- 6 Af 5 Af 

-0.42388848 - 0.35453926 Af - 0.079847022 Af 2 - 0.004109703 Af 3 + 0.00008739999 Af 4 
+2.4368878 10~ 6 Af 5 - 0.00096920612 Af Af - 0.00074246514 M 2 N - 0.00013828482 Af 3 Af 

-5.8095838 10~ 6 Af 4 Af - 5.6634686 10~ 7 Af 5 Af 

-0.057509877 - 0.048617439 Af - 0.011215648 Af 2 - 0.00061966388 Af 3 + 0.000011640607 Af 4 
+1.8658758 10- 7 Af 5 



In this way the Borel transform becomes a series in powers of y(g) that converges for all 
positive values of g provided that all singularities of the Borel transform are on the real 
negative axis []99|. In this case one obtains a convergent sequence of approximations for the 
original quantity. For the O(JV) -symmetric theory accurate estimates (see e.g. Ref. PD| ) 
have been obtained resumming the available series: the /3-function ]98[] is known up to six 
loops, while the functions 77^ and r\ t are known to seven loops [ |11U[ . 

The large-order behavior of the perturbative expansions in the MN model can be studied 
by employing the same methods used in the standard <f) A theory [|1 1 1|] . We may consider the 
series in u and v at fixed ratio z = v/u. The large-order behavior of the resulting expansion 
in powers of u is determined by the singularity of the Borel transform that is closest to the 
origin, Ub(z), given by 
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TABLE IV. The coefficients , cf. Eg. (513) 



i,j 


MN M ij 


3 , 


0.64380517 + 0.05741276 M N - 0.0017161966 Af 2 N 2 


2 , 1 


1.3928409 + 0.29248953 M + 0.0061625366 MN - 0.0030911252 M 2 N 


1 , 2 


1.0440961 + 0.26681908 M + 0.0029142164 M 2 


, 3 


0.27385517 + 0.075364029 M + 0.0018504016 Af 2 


4 , 


-0.76706177 - 0.089054667 M N + 0.000040711369 Af 2 N 2 - 0.000087586118 Af 3 N 3 


3 , 1 


-2.2398975 - 0.49868661 M - 0.043414868 MN - 0.0058040064 M a N + 0.00021024326 M z N* 
-0.00023647795 Af 3 N 2 


2 , 2 


-2.5589671 - 0.75680382 M - 0.031949535 M 2 + 0.0058629697 MN + 0.0018447227 M 2 N 
-0.00016585294 Af 3 N 


1 , 3 


-1.3553512 - 0.42919212 M - 0.022689592 M* + 0.000045447676 M 3 


, 4 


-0.27925724 - 0.091833749 M - 0.0054595646 M 2 + 0.000023722893 M 3 


5 , 


1.0965348 + 0.15791293 MN + 0.0023584631 Af 2 N 2 - 0.000061471346 Af 3 N 3 
-5.3871247 10" 6 M 4 N 4 


4 , 1 


4.0438017 + 0.94274678 M + 0.1469466 MN + 0.028781322 AP N - 0.0016180434 Af 2 N 2, 
-0.0004537935 M 3 N 2 + 7.2782691 10~ 7 M 3 N 3 - 0.000020110739 Af 4 N 3 


3 , 2 


6.2251917 + 2.0216275 M + 0.11770913 M 2 + 0.0079585974 M N - 0.002803581 M 2 N 
-0.0011929602 M 3 N + 0.00019299342 M 2 N 2 + 0.000046931931 M 3 N 2 - 0.000026294129 M 4 N 2 


2 , 3 


4.9862586 + 1.7772486 M + 0.13144724 M 2 - 0.00035323931 M A - 0.017641691 M N 
-0.0052142607 M 2 N - 0.00021963586 M 3 N - 0.000011811618 M 4 N 


1 , 4 


2.0658132 + 0.75909418 M + 0.060829135 M 2 + 0.000053192889 M A + 2.0293989 10~ o Af 4 


, 5 


0.35174477 + 0.13242502 M + 0.011322026 M 2 + 0.000054833719 M 3 + 8.6768933 10~' M 4 


ft n 
u , U 


1 r 7'7j1KCQQ H QH-IH/IQI ft M M Ci nflO A '3 Q on" 7 !'! \ 1 - 1 f\ 111111111 '(■'(')< i '- 1 v <-': 1 1 7 ■ A Af-A 

— 1.//40OOO — U.oU4U4olo 1VL iv — U.UUy4oooU / y ±VL iv + U.UUUUooyyoo04 lVl iv 

-6.5724895 10" 6 M 4 N 4 - 3.753114 10" 7 M 5 N 5 


5 , 1 


-7.9179198 - 1.9119099 M - 0.4354995 M N — 0.09835005 M 2 N + 0.0016283562 M 2 N 2 
+U.UUUO 1 4ybyUl 1V1 iv — U.UUUUoylb4oo ( M iv — U.UUUU41oUoDoO 1VL iv 
-8.0625922 1CT 7 M 4 N 4 - 1.7896837 10 -6 M 5 N 4 




1 e QKft/inc; C Q ; |'J7fi|ft 11 n Q'7'~?Q'7nftft Lf^ n 1 ftQOft AT \r D HQIOO^Vftl A f ^ /\f 

+0.00080205375 M 3 N + 0.00036992499 Af 2 Af 2 - 0.00027756993 M 3 N 2 - 0.00010718274 Af 4 A^ 2 
+8.9627582 10" 6 M 3 N 3 - 1.1514477 10" 6 A/ 4 N 3 - 3.3124581 10" 6 M 5 N 3 


3 , 3 


-16.500282 - 6.4536377 M - 0.60891787 Af 2 - 0.0068863707 Af 3 + 0.067945799 M N 
+0.025945701 M 2 N + 0.0019347624 M 3 N - 0.00010939543 Af 4 N - 0.00064953302 M 2 N 2 
-0.00017752901 M 3 N 2 - 9.1250241 10" 6 A/ 4 Af 2 - 2.8529338 10" 6 M 5 N 2 


2 , 4 


-10.296588 - 4.1925497 M - 0.43192146 M 2 - 0.0069161245 Af 3 - 0.000022535528 A/ 4 
+0.036155315 M N + 0.011884073 M 2 N + 0.00065955918 M 3 N - 0.000016702575 Af 4 N 
-9.4492956 10 -7 Af 5 N 


1 , 5 


-3.5159823 - 1.4553502 M - 0.15565998 Af 2 - 0.0028818538 M 3 + 2.4768276 10" b Af 4 
+1.2194955 10~ 7 M 5 


, 6 


-0.51049889 - 0.21485252 M - 0.023839375 A/ 2 - 0.00050021682 M 3 + 2.0167763 10"° M 4 
+4.4076733 10" 8 M 5 



1 



Ub{z) 

1 

u b {z) 



-a (Rmn + Rmz) 
J {Rmn + jjRmz 



for 
for 



z > and z < 
2N R MN 



2N R 



MN 



N+ 1 R M 



N+l R 



< z < 0. 



M 



where 



0.14777422..., 



R 



K 



+ K 



[2.22) 



(2.23) 



Using Eq. ( |2.22|) and the conformal mapping (|2.21 ), one can resum the perturbative series in 



u at fixed z. This method has been applied in Ref. [9^] to the analysis of the renormalization- 
group functions of the three-dimensional cubic model. 

The result ( gjp has been obtained for integer M } N > 1. For N = 0, one may think 
that the correct behaviour is obtained by analytic continuation of (|2.22|) , i.e. 
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TABLE V. The coefficients e\]> , cf. Eq.(^). 



i,j 


MN M c ij 


3 , 


0.00054176134 + 0.00033860084 MN + 0.000033860084 M 2 N 2 


2 , 1 


0.001625284 + 0.00081264201 M + 0.0002031605 MN + 0.00010158025 M 2 JV 


1 , 2 


0.001625284 + 0.0010158025 M + 0.00010158025 Af 2 


, 3 


0.00054176134 + 0.00033860084 M + 0.000033860084 M 2 


4 , 


0.00099254838 + 0.00070251807 M JV + 0.0001018116 M 2 N 2 - 6.5516886 10" ' M 3 N 3 


3 , 1 


0.0039701935 + 0.0019850968 M + 0.00082497551 MN + 0.00041248776 M 2 N - 5.2413509 10"° M 2 N 2 
-2.6206755 10" 6 A/ 3 JV 2 


2 , 2 


0.0059552903 + 0.0039994684 M + 0.00051091164 M 2 + 0.00021563998 MN + 0.000099957964 M 2 N 
-3.9310132 10" 6 M 3 JV 


1 , 3 


0.0039701935 + 0.0028100723 M + 0.0004072464 AI 2 - 2.6206755 10~ b M 3 


, 4 


0.00099254838 + 0.00070251807 M + 0.0001018116 M 2 - 6.5516886 10"' M 3 


5 , 


-0.00036659735 - 0.0002572117 M N - 0.000032026611 M 2 N 2 + 2.2430702 10" b M 3 N 3 
-1.1094045 10" 7 M 4 JV 4 


4 , 1 


-0.0018329867 - 0.00091649336 M - 0.00036956513 MN- 0.00018478256 M 2 N + 0.000024649511 M 2 N 2 
+0.000012324756 A/ 3 JV 2 - 1.1094045 10" 6 A/ 3 JV 3 - 5.5470225 10" 7 M 4 N 3 


3 , 2 


-0.0036659735 - 0.0024800846 M - 0.00032354892 M 2 - 0.000092032422 MN - 7.9092305 10~' M 2 N 
+0.000022612644 M 3 N + 4.0737345 10" b M 2 N 2 - 1.8194176 10" 7 A/ 3 N 2 - 1.1094045 10" b M 4 N 2 


2 , 3 


-0.0036659735 - 0.0026040248 M - 0.00035159322 M 2 + 0.000016962916 M 3 + 0.000031907838 MN 
+0.00003132711 M 2 N + 5.4677866 10~ 6 M 3 N - 1.1094045 10 _b M 4 N 


1 , 4 


-0.0018329867 - 0.0012860585 M - 0.00016013305 M 2 + 0.000011215351 M 3 - 5.5470225 10" ' Af 4 


, 5 


-0.00036659735 - 0.0002572117 M - 0.000032026611 M 2 + 2.2430702 10~ b M 3 - 1.1094045 10"' AI 4 


, U 


u.uuuoyoDoUo / -f- u.uuuoDOooy4i ivi is u.uuuizuo ( ouz ivi iv o. 1 4Doy /y iu ivi i\ 
-3.8385183 10~ 8 A/ 4 N 4 - 1.0441273 10" s M 5 JV 5 


5 , 1 


0.0041740822 + 0.0020870411 M + 0.0013081154 MN + 0.00065405768 M 2 N + 0.000069380438 M 2 N 2 
+U.UUUUo4uyUziy Ixi Is — Z.lUUol04 1U LVI IV — l.UuUlOoz 1U IVI IV — i.zozyozo 1U IVI IV 
-6.264764 10~ 8 Af 5 JV 4 




U.UlU4oozUo H~ U.UU ( 4ZU4ozy ikz -j- U.UU11U144U1 iVi + U.UU1UD /4Uoo ikz iV + U.UUUoyo4yoo4 iw iv 

+0.00008239735 M 3 JV + 8.6563947 10" 6 M 2 N 2 + 3.9934653 10" 6 M 3 N 2 - 1.6736603 10" 7 M 4 N 2 
-1.9034703 10" 7 M 3 N 3 - 4.0841171 10" 7 M 4 N 3 - 1.566191 10" 7 M 5 N 3 


3 , 3 


0.013913607 + 0.010984227 M + 0.0021785508 AP + 0.000082419651 A/ 3 + 0.00033296155 M N 
+0.00023139203 M 2 N + 0.000031719215 M 3 N - 3.6820589 10 -7 AI 4 N + 1.5175718 10 -6 M 2 N 2 
+7.9509223 10" 7 M 3 N 2 - 3.9949777 10" 7 Af 4 TV 2 - 2.0882547 10" 7 A/ 5 /V 2 


2 , 4 


0.010435206 + 0.0084674045 M + 0.001790972 AP + 0.000082410852 M 3 - 3.1236483 10" ' M 4 
+0.000020486706 M N + 0.000017623286 Af 2 JV + 3.7896168 10~ 6 AI 3 JV - 2.6341291 10" 7 M 4 N 
-1.566191 10 _7 M 5 JV 


1 , 5 


0.0041740822 + 0.0033951565 M + 0.00072343812 M 2 + 0.000034480187 AI 3 - 2.303111 10" ' M 4 
-6.264764 10~ 8 M 5 


, 6 


0.00069568037 + 0.00056585941 M + 0.00012057302 AI^ + 5.7466979 10"° M 3 - 3.8385183 10" 5 M 4 
-1.0441273 10- 8 M 5 



u b lz 



-a (| + R M z) 



(2.24) 



for all z. However, this is not correct. Indeed, as explicitly shown in Refs. ||101|J102|| in the 
context of the zero-dimensional random Ising model, there is an additional contribution to 
the large-order behavior of the series in u at fixed A = —u/v, which makes the series non- 
Borel summable, giving rise to singularities of the Borel transform on the positive real axis. 
They are due to the zeroes of the partition function at fixed disorder. We have no reason to 
believe that similar non-Borel summable contributions are not present in higher dimensions. 
It is likely that the same phenomenon occurs even in three dimensions. As a consequence, 
a summation procedure based on Eq. (|2.24 ) and a conformal mapping of the type Q2.21 
would not lead to a sequence of approximations converging to the correct result [|101 
Fortunately, this is not the end of the story. As shown recently in Ref. 



103| 1, at least 



in zero dimensions, one can still resum the perturbative series. Indeed the zero- dimensional 
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TABLE VI. The coefficients eg j , cf. Eg. (gig). 



i 




MN M c ij 


3 





-0.025120499 - 0.016979919 M N — 0.0022098349 M z N z 


2 


i 


-0.075361497 - 0.037680748 M - 0.013259009 MN- 0.0066295047 M z N 


1 


2 


-0.075361497 - 0.049233615 M - 0.0057764331 M z - 0.0017061432 MN - 0.00085307161 M 2 N 





3 


-0.025120499 - 0.016979919 M - 0.0022098349 M z 


4 





0.021460047 + 0.015690833 M N + 0.0024059273 M 2 N 2 - 0.000037238563 M 3 N 3 


3 


1 


0.08584019 + 0.042920095 M + 0.019843236 M N + 0.0099216178 M 2 N - 0.0002979085 M z N 2 
-0.00014895425 M 3 N 2 


2 


2 


0.12876028 + 0.087481041 M + 0.011550449 M z + 0.0066639546 M N + 0.0028851146 M' z N 
-0.00022343138 M 3 N 


1 


3 


0.08584019 + 0.061713883 M + 0.0089678045 M' z - 0.00021454473 M 3 + 0.0010494477 M N 
+0.00065590478 M 2 N + 0.000065590478 M 3 N 





4 


0.021460047 + 0.015690833 M + 0.0024059273 M' z - 0.000037238563 M 3 


5 





-0.022694287 — 0.017985168 M N — 0.0035835384 M z N z - 0.00013566164 M 3 N 3 
-1.699309 10" 6 M 4 N A 


4 


1 


-0.11347143 - 0.056735717 M - 0.033190124 M N - 0.016595062 M' z N - 0.0013226302 M' z N 2 
-0.00066131512 M 3 N 2 - 0.00001699309 M 3 N 3 - 8.496545 10" 6 M 4 N 3 


3 


2 


-0.22694287 - 0. 15931239 M - 0.022920477 M z - 0.020539294 MN - 0.012755939 M' z N 
-0.001243146 M 3 N - 0.00015896839 M 2 N 2 - 0.00011347037 M 3 N 2 - 0.00001699309 M 4 N 2 


2 


3 


-0.22694287 - 0.1746255 M - 0.032245328 M z - 0.00083414751 M 3 - 0.0052261822 M N 
-0.0035900565 M 2 N - 0.00052246891 M 3 N - 0.00001699309 M 4 N 


1 


4 


-0.11347143 - 0.089150636 M - 0.01739473 M' z - 0.00062322789 M 3 - 0.000014796237 M 4 
-0.0007752054 MN - 0.0005229621 M 2 N - 0.000055080316 M 3 N + 6.2996919 10~ 6 M 4 N 





5 


-0.022694287 - 0.017985168 M - 0.0035835384 M z - 0.00013566164 M 3 - 1.699309 10~ b Af 4 





U 


n non/ic;nfiin i n no i 07 1 trvn a t at i n nnt^voQQO'? v- 1 n nnnQ 1 kk^qrq r r3 at<$ 

u.uzy4ouoiy + u.uz4o / 4o / y ivi i\ u.uuo < zooy i ivi i\ -t- u.uuuoioo (ooo ivi i\ 
-5.858689 10~ 6 M 4 TV 4 - 1.0373506 10~ 7 M 5 TV 5 


5 


1 


0.17670371 + 0.088351856 M + 0.060895618 M N + 0.030447809 M z N + 0.0039225729 M z N' z 
+U.UUiyolZou4 M iv — U.UUUUo 1 ol4oz / IVL iv — U.UUUUooyU / olo M iV — l.z44ozUo 1U ili iv 
-6.2241039 10" 7 M 5 TV 4 


4 


2 


0.44175928 + 0.31708477 M + 0.048102565 M' z + 0.056033914 M N + 0.037304669 M' z N 
-(-u.uu4d4ooooo ivi iv + u.uuuoio / iyoo ivi iv + u.uuuuyo^u / do / iw i\ — u.uuuuozudiud ivi ±\ 
-5.434447 10~ 6 A/ 3 - 5.8292754 10~ 6 A/ 4 W 3 - 1.556026 10~ 6 M 5 N 3 


3 


3 


0.58901238 + 0.47310641 M + 0.097561547 M 2 + 0.0041307177 M 3 + 0.024385168 M N 
+0.017002488 M 2 N + 0.002190569 M 3 N - 0.00010719141 M 4 N + 3.9039471 10" 6 M 2 N 2 
-9.7139715 10" 6 M 3 N 2 - 9.9823751 10" 6 Af 4 N 2 - 2.0747013 10~ 6 M 5 N 2 


2 


4 


0.44175928 + 0.36731802 M + 0.081706647 M z + 0.0041012326 M 3 - 0.000071248084 M 4 
+0.0058006649 M N + 0.0042193071 M 2 N + 0.00063244692 M 3 N — 0.000016632252 M 4 N 
— 1.556026 10 -6 M 5 N 


1 


5 


0.17670371 + 0.14827827 M + 0.033627917 M 2 + 0.001755187 M 3 - 0.000040961717 M 4 
-1.1887572 10" 6 M 5 + 0.00096920611 MN+ 0.00074246513 M 2 N + 0.00013828482 M 3 N 
+5.809583 10" 6 Af 4 N + 5.6634686 10~ 7 M 5 N 





6 


0.029450619 + 0.024874579 M + 0.005728397 M 2 + 0.00031557863 M 3 - 5.858689 lO" M 4 
-1.0373506 10~ 7 M 5 



free energy can be obtained from its perturbative expansion if one applies a more elaborated 
procedure which is still based on a Borel summation. Let us write the double expansion of 
the free energy f(u, v) in powers of u, v as 



f(u,v) = J2 c «( 



v 111 



n=0 

oo 



k=0 



(2.25) 
(2.26) 



The main result of Ref. ||103|| is that the expansions of the coefficients ( |2.26|) and the resulting 
series at fixed v, Eq. ( 2.25|) , are Borel summable. Using this result, a resummation of the 
free energy is obtained in two steps. First, one resums the coefficient c n (v); then, using the 
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computed coefficients, one resums the series in u. The resummation of Eq. ( 2.26|) can be 
performed using the Pade-Borel-Leroy method, as suggested in Ref. ||103|| . However, also 
the conformal method can be used, since the large-order behavior is known exactly. Indeed, 



[v) oc 



d n f(u, v) 



[2.27) 



u=0 



Thus, c n (v) can be related to zero-momentum correlation functions in the theory with u = 0, 
which is the standard M-vector model. Therefore, one can use the well-known results for 
the large-order behavior of the perturbative series in the 0(M)-symmetric theory [|9~6"|,|97"|. 



III. ANALYSIS OF THE SIX-LOOP EXPANSION FOR N = 



A. The random Ising model 



As we said in the Introduction, the random Ising model corresponds to M — 1 and 
N = 0. There are two relevant fixed points, the Ising and the random point, see Fig. |l|. 
In Ref. [92] we already discussed the stability of the Ising point. We found that this fixed 
point is unstable since the stability matrix has a negative eigenvalue u = —0.177(6), in 
good agreement with the general argument predicting uj = —OLijvi = —0.1745(12). We will 
now investigate the random fixed point, which is stable and determines therefore the critical 
behaviour of the RIM. 

In order to study the critical properties of the random fixed point, we used several differ- 
ent resummation procedures, according to the discussion of the previous Section. Following 
Ref. [ |103]| , for each quantity we consider, we must perform first a resummation of the series 
in v, see Eq. (|2.26|) . This may be done in two different ways. We can either use the Pade- 
Borel method, or perform a conformal mapping of the Borel-transformed series, using the 
known value of the singularity of the Borel transform. Explicitly, let us consider a p-loop 
series in u and v of the form 



p p—n 



E E CnkU 

n=0 fc=0 



V. 



(3.1) 



In the first method, for each < n < p, we choose a real number b n and a positive integer 
r n such that r n < p — n; then, we consider 



Ri(c n )(p; b n ,r n ;v) 



dt e~H 



tj> n 



'p-n-r n 

E 

. i=0 



i+Ea(ft 



i=l 



(3.2) 



The coefficients B { and Cj are fixed so that Ri(c n )(p; b n , , ,,. 



n-= cn k v k + 0(vv- 



Here we are resumming the Borel transform of each coefficient of the series in u by means 
of a Pade approximant [p — n — r n /r n ]. Eq. ( |3.2|) is well defined as long as the integrand 
is regular for all positive values of t. However, for some values of the parameters, the Pade 
approximant has poles on the positive real axis — we will call these cases defective — so 
that the integral does not exist. These values of b n and r n must of course be discarded. 
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The second method uses the large-order behaviour of the series and a conformal mapping 
99[Jl00fl. In this case, for each < n < p, we choose two real numbers b n and a n and consider 



R2{c n ){p; b n , r n , v)=Y,B k dte-H K —^ 



(y(vt)Y 



y(vt)) a "' 



(3.3) 



where 



y(t) 



9i\+t 



9i\+t + 



(3.4) 



and Ijj = — 1/a is the singularity of the Borel transform for the pure Ising model. The 
numerical value of a is given in Eq. ( 2.23|) . Using these two methods we obtain two different 
partial resummations of the original series ( |3.1| ): 



Ri(c n )(p; b n ,r n ;v)u n , 



n=0 
V 



E R 2(cn) (p; b n , a n ; v)u n . 

n=0 



(3.5) 
(3.6) 



Nothing is known on the asymptotic behaviour of these series, and we will thus use the 
Pade-Borel method. Starting from Eq. Q3.5D we will thus consider 



poo 

Ei(c)(g,p;6 n ,r u ;{6 n },{r n }) = / dte'H K 

Jo 



q-r u 

E 

. i=0 



E B t (v)(t u y 



i + J2c t (v)(t u y 



i=l 



(3.7) 



The coefficients Bi{v) and Ci{v) are fixed so that Ei(c)(q,p; b u , r u \ {b n }, {r n }) coincides with 
the expansion (|3.5|) up to terms of order 0(u q+l ). Note that we have introduced here three 
additional parameters: b u , the power appearing in the Borel transform, r u that fixes the order 
of the Pade approximant, and q which indicates the number of terms that are resummed, and 
that, in the following, will always satisfy q < p — 1. Analogously, starting from Eq. ( |3.6| ), we 
define E2(c)(q,p; b u ,r u ; {b n }, {a n }). We will call the first method the "double Pade-Borel" 
method, while the second will be named the "conformal Pade-Borel" method. 

Let us now apply these methods to the computation of the fixed point (u*,v*). In 
this case, we resum the /3-functions (3 u /u and (3 v /v and then look for a common zero with 
u < 0. We consider first the resummation Ei, Eq. (|3.7| ). A detailed analysis shows that the 
coefficients R\ can only be defined for r n — 1. We have also tried r = 2 and r\ = 2, but 
the resulting Pade approximants turned out defective. Therefore, we have fixed r n = 1 for 
all < n < 5. We must also fix the parameters {&„}. It is impossible to vary all of them 
independently, since there are too many combinations. For this reason, we have taken all b n 
equal, i.e. we have set b n = b v for all n. Finally, we have only considered the case q = p — 1. 
Therefore, the analysis is based on the approximants 



Ei{-)(p; b u ,r u ;b v ) = E^ip- l,p; b u ,r u ; {b n = b v }, {r n = 1}). 



(3.8) 



Estimates of the fixed point (u*,v*) have been obtained by solving the equations 
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Ei{(3u/u){p;b u ,r u ;b v ) = 0, E 1 (p 1J /v)(p; b u , r u ; b v ) = 0. (3.9) 

We have used r u = 1, 2, p = 4, 5, 6 and we have varied b u and 6^ between and 20. As 
usual in these procedures, we must determine the optimal values of the parameters b u and 
b v . This is usually accomplished by looking for values of b u and b v such that the estimates 
are essentially independent of the order p of the series. In the present case, we have not been 
able to find any such pair. Indeed, the five-loop results (p = 5) are systematically higher 
than those obtained with p = 4 and p — 6. For instance, if we average all estimates with 

< b u ,b v < 5 we obtain 

u* = -0.66(1) (p = 4), -0.78(2) (p = 5), -0.63(3) (p = 6); (3.10) 

r = 2.235(3) (p = 4), 2.273(4) (p = 5), 2.250(23) (p = 6). (3.11) 

The uncertainties quoted here are the standard deviations of the estimates in the quoted 
interval and show that the dependence on b u , b v , and r u is very small compared to the 
variation of the results with p. Increasing b u ,b v does not help, since the five-loop result is 
largely insensitive to variations of the parameters, while for p = 4 and p = 6 \u*\ and v* 
decrease with increasing b u and b v . It is difficult to obtain a final estimate from these results. 
We quote 

u* = -0.68(10), v* = 2.25(2), (3.12) 

that includes all estimates reported above. 

The instability of the results reported above with p seems to indicate that some of the 
hypotheses underlying the choice of the parameters is probably incorrect. One may suspect 
that choosing all b n equal does not allow a correct resummation of the coefficients, and 
that f3u and (3^ need different choices of the parameters. We have therefore tried a second 
strategy. First, for each /3-function, we have carefully analyzed each coefficient of the series 
in u, trying to find an optimal value of the parameter b n — r n was fixed in all cases equal to 

1 — by requiring the stability of the estimates of the coefficient with respect to a change of 
the order of the series. However, only for the first two coefficients we were able to identify a 
stable region, so that we could not apply this method. On the other hand, as we shall see, 
this method works very well for the resummations of the coefficients that use the conformal 
mapping. 

Let us now discuss the conformal Pade-Borel method. As before, we have tried two 
different strategies. In the first case we have set all b n equal to b v and all a n equal to a v , 
we have used the same parameters for the two /5-functions, and we have looked for optimal 
values of b u , b v and a v , setting r u = 1, 2. While before, for each p, the estimates were stable, 
in this case the fluctuations for each fixed p are very large, and no estimate can be obtained. 

Then, we applied the second strategy. We analyzed carefully each coefficient of the 
series in u, finding optimal values fo ni0p t and a nj0pt for each n and /3-function. Of course, the 
required stability analysis can only be performed if the series is long enough, and thus we 
have always taken q < 4. Therefore, we consider 

i P'i Til] 8b, $a + 5b}, {«n,opt + 8a}), (3.13) 
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TABLE VII. Estimates of (u*,v*) obtained using the conformal Pade-Borel method. The 
results are averages over < b u < 20, —3 < 6b < 3, — 1 < S a < 1, r 2)U = 1,2. "Def" is the 
percentage of defective Pade approximants in each analysis. 



q 


ri )U 


p = 5 

(«*, V*) 


Def 


(u*, 


p = 6 

V*) 


Def 


3 





[-0.615(27), 2.175(43)] 


80% 


[-0.641(7), 


2.199(10)] 


78% 


4 





[-0.618(22), 2.190(35)] 


74% 


[-0.630(9), 


2.194(16)] 


61% 


4 


3 


[-0.619(20), 2.191(33)] 


71% 


[-0.632(8), 


2.196(15)] 


61% 



where 5 a and 5 b are (n-independent) numbers which allow us to vary b n and a n around the 
optimal values. Estimates of the fixed point are obtained from 

E 2 (Pu/u)(q,p;b u ,r ltU ;5 b ,5 a ) = 0, E 2 ((3jj/v)(q,p;b u ,r 2:U ; 5 b , 5 a ) = 0. (3.14) 

The first problem which must be addressed is the value of the parameters r l u and r 2jU . For 
ftu/u we find that the Pade approximants are always defective for r l n = 1,2; they are well 
behaved only for r l n = 3 and q = 4. Since the resummed series in u has coefficients that 
are quite small, we decided to use also r l u = 0, which corresponds to a direct summation 
of the series in u, without any Pade-Borel transformation. For fly/v we did not observe a 
regular pattern for the defective Pade's and we have used r 2u = 1,2, discarding all defective 



cases. The results, for chosen values of p, q, and ri M are reported in Table [VII] . The quoted 
uncertainty is the standard deviation of the results when —3 < 5 b < 3 and — 1 < 5 a < 1. 
This choice is completely arbitrary, but in similar analyses of different models we found 
that varying a by ±1 and b by ±3 provides a reasonable estimate of the error. Notice that 
we have not optimized b u , but we have averaged over all values between and 20, since 
the dependence on this parameter is extremely small. The results are stable, giving a final 
estimate (average of the results with p = 6, q = 4) 

u* = -0.631(16), v* = 2.195(20). (3.15) 

The error bars have been chosen in such a way to include all central values for p = 5 and 
p = 6. It should be noted that, even if our results are quite stable with respect to changes 
of the parameters, most of the approximants do not contribute since they are defective. For 
these reasons, in the following we will always carefully check the dependence of the estimates 
on the value of the fixed point, considering also values of (u*,v*) that are well outside the 
confidence intervals of Eq. (|3.15| ). 



We can compare our results for the fixed point with previous determinations. Ref. f24 
reports (u*,v*) = (—0.667,2.244) obtained from the Chisholm-Borel analysis of the four- 
loop series. The same expansion was also analyzed by Varnashev [|35| obtaining (u*,v*) = 
(-0.582(85), 2.230(83)) and (7i*,7J*) = (-0.625(60), 2.187(56)) using different sets of Pade 
approximants. The e-algorithm by Wynn with a Mittag-Leffler transform was used in Ref. 
p5| finding (u*,v*) = (—0.587, 2.178). From the analysis of the five-loop series Pakhnin and 
Sokolov obtain (u*,v*) = (-0.711(12), 2.261(18)). While the four-loop results are in 
good agreement with our estimates, the five-loop estimate differs significantly, a fact that 



may indicate that the claim of Ref. 37 that the error on their estimates is approximately 
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1-2% is rather optimistic. Note also that the five-loop result is quite different from the 
previous four-loop estimates. 

We have also tried to determine the eigenvalues of the stability matrix Q, cf. Eq. ( |2 . 8| ) , 
that controls the subleading corrections in the model. We used both a double Pade-Borel 
transformation and the conformal-Pade-Borel method. In the first case we obtain estimates 
that vary strongly with the order, and, as it happened for the position of the fixed point, it 
is impossible to obtain results that are insensitive to the order p. For r u — 1, discarding the 
cases in which the computed eigenvalues are complex, we obtain for the smallest eigenvalue 
to: 

u = 0.16(2) (p = 4), 0.21(3) (p = 5), 0.16(3) (p = 6). (3.16) 

We have included in the error the dependence on the position of the fixed point. These 
estimates have been obtained setting r n = 1 and averaging over b u and b v varying between 
and 10. We have not tried to optimize the choice of these parameters, since the estimates 
show only a small dependence on them. We have also considered r u = 2. In this case a large 
fraction of the approximants is defective (for p = 4 they are all defective). We obtain 

u = 0.29(4) (p = 5), 0.33(5) (p = 6). (3.17) 

The quite large discrepancy between the estimates Q3.16]) and ( |3.17| ) clearly indicates that 
the analysis is not very robust. A conservative final estimate is 

oj = 0.25(10), (3.18) 

that includes the previous results. 

We have also tried the conformal-Pade-Borel method, optimizing separately each coef- 
ficient . However, several problem appeared immediately. First, we could not perform a 
Pade-Borel resummation of the series in u of the elements of the stability matrix. Indeed, in 
all cases, some Pade approximant was defective. As in the determination of the fixed-point 
position, we tried to resum the series in u without any additional transformation. For p = 4 
this gives reasonable results, and we can estimate u = 0.29(9). However, for p = 5, 6 all 
eigenvalues we find are complex, and as such must be discarded. 

The fact that the series appearing in the stability matrix generate always defective Pade 
approximants may indicate that the series in u are not Borel summable. In this case, 
one expects that the estimates converge towards the correct value up to a certain number 
of loops. Increasing further the length of the series, worsens the final results. If indeed 
the expansion is not Borel summable, the previous results seem to indicate that for the 
subleading exponent oj the best results are obtained at four loops. 

Let us now compute the critical exponents. As before, we tried several different methods. 
A first estimate was determined using the double Pade-Borel method. Each exponent was 
computed from the approximants i?i(e)(p; b u , r u ; b v ) defined in Eq. (|3.8|). For 7 and u, the 
series I/7 and 1/v are more stable and thus the final estimates are obtained from their 
analysis. For rj, if we write rj = Vn{v)u n , then rjo(v) ~ v 2 and r)i(v) ~ v. In this case, for 
n = we resummed the series t]q(v)/v 2 , while for n = 1 we considered r)i(v)/v. The results 
we obtain are very stable, even if we do not optimize the parameters b u and b v . Without 
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TABLE VIII. Estimates of the critical exponents for the RIM universality class obtained using 
a double Pade-Borel resummation. The estimates correspond to the following choices of the pa- 
rameters: < b u , b v < 20, r u = 1, 2. The subscript indicates the percentage of cases in which some 
Pade approximant appearing in the resummation procedure was defective. 



(U* v*\ 


7 i 






7 


'/ 






n — 
P — 


4 






( n fi^i 9 i 

I — U.UOl, £.laOJ 

(-0.640,2.230) 


u.u / uu^o j 49% 
0.6740(3) 49% 






1.32 82(5) 50% 




(-0.680,2.240) 


0.6725(2) 49% 






1.3283(5)50% 








P = 


5 






(-0.631,2.195) 


0.6740(4) 0% 






1.3270(7)o% 


0.03081(3) % 


(-0.745,2.275) 


0.6739(4) 0% 






1.3291(7)o% 


0.02738(2) 0% 


(-0.800,2.270) 


0.6687(3) 0% 






1.3213(7)o% 


0.02387(2) 0% 






P = 


6 






(-0.631,2.195) 


0.6677(4) % 






1.3130(8)o% 


0.03272(2) 50 % 


(-0.570,2.210) 


0.6745(6) 0% 






1.3225(11) % 


0.03810(4) 50% 


(-0.700,2.290) 


0.6727(4) 0% 






1.3223(8)o% 


0.03333(2) 50 % 



choosing any particular value for them, but simply averaging over all values between and 
10, we obtain the results of Table |VTII| . Note that we have not quoted any estimate of r] 
for p = 4: in all cases, some Pade approximant was defective. The quoted uncertainty, that 
expresses the variation of the estimates when changing b u , b v , and r u , is very small, and it is 
clear that it cannot be interpreted as a correct estimate of the error, since the variation with 
the order p of the series is much larger. In Table |V11J| we also report the estimates of the 
exponents corresponding to several different values of (u*,v*): beside the estimate ( J3.15|) , we 
consider two values appearing in the first analysis of the fixed point position, those with the 
largest and smallest value of u*, when b u and b v vary in [0,5]. The dependence on (u*,v*) 
is quite small, of the same order of the dependence on the order p. As final estimate we 
quote the value obtained for p = 6, using the estimate ( 3.15Q for the fixed point. The error 



is estimated by the difference between the results with p = 5 and p = 6. Therefore we have 
7 = 1.313(14), v = 0.668(6), r? = 0.0327(19). (3.19) 



Note that, within one error bar, all estimates of 7 and v reported in Table [VII]] are compatible 



with the results given above. Instead, the estimates of r\ show a stronger dependence on 
the critical point, and a priori, since we do not know how reliable are the uncertainties 
reported in Eq. ( [3.15 ), it is possible that the correct estimate is outside the confidence 



interval reported above. 

We will now use the conformal Pade-Borel method. A first estimate is obtained consid- 
ering approximants of the form 

E 3 (-)(q,p;b u ,r u ;b v ,a v ) = E 2 (-)(q,p;b u ,r u ; {b n = b v },{a n = a v }), (3.20) 

setting all b n equal to b v and a n equal to a v . The results show a tiny dependence on 
b u , while no systematic difference is observed between the approximants with r u = 1 and 
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TABLE IX. Estimates of the critical exponents for the RIM universality class obtained using 
the conformal Pade-Borel method (first analysis). All estimates correspond to the following choices 
of the parameters: < b u < 10, 2 < b v < 8, —1.5 < a v < 0.5, r u = 1,2. The subscript indicates 
the percentage of defective Pade approximants. The upper part of the table reports the estimates 
for (u*,v*) = (—0.631,2.195) and different values of p and q. The lower part reports estimates 
obtained setting p = 6 and q = 4 for several values of (u*,v*). 





p = 5, q = 3 


p = 5, q = 4 


p = 6, q = 3 


p = 6, q = 4 


7 


1.336(18)40% 


1.308(19)40% 


1.305(15) 50 % 


1.338(21)o% 


V 


0.673(11)44% 


0.660(10) 44% 


0.657(7) 54 % 


0.676(ll)o% 


V 


0.0285(15) 4% 


0.0288(17) 46% 


0.0278(4) 27 % 


0.0279(5) 31 % 




(-0.647,2.215) 


(-0.615,2.175) 


(-0.700,2.290) 


(-0.570,2.210) 


7 


1.342(21) 0% 


1.334(21)o% 


1.371(23)o% 


1.349(21)o% 


V 


0.678(11) 0% 


0.674(11) 0% 


0.694(13)o% 


0.685(13)o% 


V 


0.0277(5) 31% 


0.0280(4) 33% 


0.0270(12) 28% 


0.0334(4) 30% 



r u = 2. Therefore, we averaged over all non-defective results with < b u < 10 and r u = 1,2. 
Then, we looked for optimal intervals [a opt — Aa, a opt + Aa], [6 opt — Ab, b opt + Ab] for the 
parameters a v and b v . They were determined by minimizing the discrepancies among the 
estimates corresponding to (p, q) = (5, 3), (5, 4), (6, 3), and (6, 4). Using Aa = 1 and Ab = 3 
as we did before, we obtain 6 opt = 5 and a opt = —0.5. The results corresponding to this 
choice of parameters are reported in Table [IX|. As final estimate we quote the value obtained 
for p = 6 and q = 4, using the estimate ( |3.15|) for the fixed point: 

7 = 1.338(21), v = 0.676(11), rj = 0.0279(5). (3.21) 

For 7 and v the estimates given above are compatible with all results of Table [IX|. In 
particular, they are correct even if the error in Eq. ( [3.15 ) is underestimated. They are also 



in good agreement with the estimates obtained with the double Pade-Borel transformation, 
cf. Eq. ( p,19|) . On the other hand, it is not clear if the error on i] is reliable. Indeed, 



comparison with Eq. ( |3.19| ) may indicate that the correct value of r] is larger than what 
predicted by this analysis. 

As we did for the fixed point, we can also use the approximants E 2 denned in Eq. ( p,13|) , 
optimizing separately each coefficient. The results are reported in Table [X] and correspond 
to < b u < 10, —3 < 5b < 3, — 1 < 5 a < 1 and r u = 1, 2. As it can be seen from the 
very small "errors" on the results, the dependence on b u is tiny and we have not tried to 
optimize this parameter. The results are reasonably stable with respect to changes of p and 
q and also the dependence on the value of the fixed point is small. As final results from this 
analysis we quote the values obtained with p = 6 and q = 4: 

7 = 1.321(8), v = 0.681(1), 77 = 0.0313(5). (3.22) 



We can compare these results with the previous estimates ( |3.19| ) and ( p.21[ ). The agreement 



is reasonable, although the quoted error on v and 77 is probably underestimated. This is 
confirmed by the fact that the scaling relation 7 = (2 — 77)// is not satisfied within error bars: 
indeed, using the estimates of v and 77, we get 7 = 1.341(2). 
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TABLE X. Estimates of the critical exponents for the RIM universality class obtained using 
the conformal Pade-Borel method (second analysis). All estimates correspond to the following 
choices of the parameters: < b u < 10, — 3 < <5& < 3, — 1 < 5 a < 1, r u = 1,2. The subscript 
indicates the percentage of defective Pade approximants. The upper part of the table reports the 
estimates for (u*,v*) = (—0.631,2.195) and different values of p and q. The lower part reports 
estimates obtained setting p = 6 and q = 4 for several values of (u*,v*). 





p = 5, q = 3 


p = 5, q = 4 


p = 6, q = 3 


p = 6, q = 4 


7 


1.335(16)42% 


1.329(15) 67% 


1-322(7)47% 


1.321(8)56% 


V 


0.684(3) 45% 


0.682(3) 72% 


0.682(1) 33% 


0.681(1) 67% 


V 


0.0299(4) 0% 


0.0299(15) % 


0.0312(6) 0% 


0.0313(5) 0% 




(-0.647,2.215) 


(-0.615,2.175) 


(-0.700,2.290) 


(-0.570,2.210) 


7 


1.323(8) 67% 


1.318(8) 67% 


1.332(10) 67% 


1.331(7) 67% 


V 


0.683(1) 68% 


0.680(1) 65% 


0.690(1) 71% 


0.687(1) 68% 


V 


0.0314(5) 0% 


0.0313(5) 0% 


0.0315(5) 0% 


0.0368(6) 0% 



We want now to obtain final estimates from the analyses given above. Since in the 
conformal Pade-Borel method we make use of some additional information, the position of 
the singularity of the Borel transform, we believe this analysis to be the more reliable one. 
As our finale estimate we have therefore considered the average between ( |3.21| ) and ( |3.22| ), 



fixing the error in such a way to include also the estimates (|3.19|) . In this way we obtain 



7 = 1.330(17), v = 0.678(10), rj = 0.030(3). (3.23) 

A check of these results is provided by the scaling relation 7 = 1/(2 — 1]). Using the values of v 
and 77 we obtain 7 = 1.336(20) in good agreement with the direct estimate. A second check 
of these results is given by the inequalities v > 2/3 ~ 0.66667 and 7 + 277/3 > 4/3 ~ 1.33333, 
that are clearly satisfied by our estimates, e.g. 7 + 2r//3 = 1.350(17). Finally, we want to 
stress that our final estimates ( |3.23| ) are compatible with all results appearing in Tables 
VIII| , p(| , and |X|, even those computed for (u*, v*) largely different from the estimate fl3.15|) . 



Thus, we believe that our error estimates take properly into account the uncertainty on the 
position of the fixed point. 

Using the scaling relations a = 2 — 3z/ and (3 = \v{\ + 77) we have 

a = -0.034(30), (3 = 0.349(5). (3.24) 

For comparison we have also performed the direct analysis of the perturbative series, 
resumming the expansions for fixed v*/u*. In zero dimensions, these series are not Borel 
summable, and this is expected to be true in any dimension. However, for the short series 
we are considering, we can still hope to obtain reasonable results. We have used the same 
procedures described in Ref. ||92|| , performing a conformal transformation and using Ub(z) 
given in Eq. ( |2.24j ) as position of the singularity. We obtain 

u* = -0.763(25), v* = 2.306(43); (3.25) 

7 = 1.327(12), v = 0.673(8), 77 = 0.029(3), u = 0.34(11). (3.26) 
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TABLE XI. Field-theoretic estimates of the critical exponents for the RIM universality class. 
Here "d = 3 exp." denotes the massive scheme in three dimensions, "MS exp." the minimal 
subtraction scheme without e-expansion. All perturbative results have been obtained by means 
of Pade-Borel or Chisholm-Borel resummations, except the results of Ref. [25] indicated by "eW" 
obtained using the e-algorithm of Wynn and of this work. 





Method 


7 


V 


V 


uJ 


This work 


d = 3 exp. 0(g 6 ) 


1.330(17) 


0.678(10) 


0.030(3) 


0.25(10) 



Ref. 
Ref. 

Ref. 



2000 
2000 

1999 



d 
d 



3 exp. 0{g 5 ) 
3 exp. 0(g 4 ) 



Ref. p§, 1989 



Ref. H, 1989 
Ref. pl, 1989 



d = 3 exp. 0(g 4 ) 
d = 3 exp. 0(g 4 ) 
d = 3 exp. 0(g 4 ) eW 
d = 3 exp. 0(g 4 ) 



1.325(3) 
1.336(2) 
1.323(5) 

1.321 
1.318 
1.326 



0.671(5) 
0.681(12) 
0.672(4) 

0.671 
0.668 
0.670 



0.025(10) 
0.040(11) 
0.034(10) 



0.034 



0.32(6) 
0.31 
0.33 

0.372(5) 



Ref. | 


,4p6), 1999 


d = 3 MS 0{g 4 ) 


1.318 


0.675 


0.049 


0.39(4) 


Ref. 


0, 1982 


scaling field 


1.38 


0.70 




0.42 



The estimate of the fixed point is very different from that computed before. This may 
indicate that the non-Borel summability causes a large systematic error in this type of 
analysis. Probably the optimal truncation for the /3-functions corresponds to shorter series. 
On the other hand, the critical exponents show a tiny dependence on the position of the fixed 
point. The estimates we obtain are in good agreement with our previous ones, indicating 
that the exponent series are much better behaved. 

Let us now compare our results with previous field-theoretic determinations, see Table 
XII. We observe a very good agrement with all the reported results. Note that our error bars 
on 7 and v are larger than those previously quoted. We believe our uncertainties to be more 
realistic. Indeed, we have often found in this work, that Pade-Borel estimates are insensitive 
to the parameters used in the analysis, in particular to the parameter b characterizing the 
Borel-Leroy transform. Therefore, error estimates based on this criterion may underestimate 
the uncertainty of the results. The perturbative results reported in Table |XJ correspond to 
the massive scheme in fixed dimension d = 3 and to the minimal subtraction scheme without 
e-expansion. It should be noted that the latter scheme does not provide any estimate at five 
loops |i6 . Indeed, at this order the resummed /3-function do not have any zero in the region 



u < 0. This fact is probably related to the fact that the series which is analyzed is not Borel- 
summable. Therefore, perturbative expansions should have an optimal truncation beyond 
which the quality of the results worsens. For the /3-functions in the minimal subtraction 
scheme, the optimal number of loops appears to be four. Two other methods have been 
used to compute the critical exponents: the scaling-field method |L7j and the A/e-expansion 
|nj . The former gives reasonable results, while the latter is unable to provide quantitative 



estimates of critical quantities, see e.g. Ref. |36]. We can also compare our results with the 



recent Monte Carlo estimates of Ref. [5~I]. The agreement is quite good for 7 and v, while 
our estimate of i] is slightly smaller, although still compatible within one error bar. This is 
not unexpected and appears as a general feature of the d = 3-expansion: indeed, also for 
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TABLE XII. Estimates of the subleading exponent ui at the 0(M)-symmetric fixed point. The 
last column reports the theoretical prediction — ccm/vm- 



M 
2~ 



p = 4 



p = 5 



p = 6 



final 



-olm/vm 



0.009(36) 



0.003(9) 



0.007(4) 



0.007(8) 0.0192(6) Ref. |l| 

0.0223(24) Ref. |8| 
0.0163(67) Ref. I§ 



3 

T 



0.142(19) 



0.151(8) 



0.156(5) 



0.156(10) 



0.172(14) Ref. M 



0.203(12) Ref. [112 



0.274(35) 



0.269(10) 



0.280(6) 



0.280(12) 



0.301(22) Ref. @ 



0.360(16) Ref. fllTl 



0.580(40) 



0.563(13) 



0.586(8) 



0.586(16) 



0.683(11) Ref. M 



the pure model, the estimate of rj obtained in the fixed- dimension expansion is lower than 
the Monte Carlo and high-temperature results (see Ref. |7(J and references therein). 



B. The random M- vector model for M > 2 



In this Section we consider the random vector model for M > 2. First, we have studied 
the region u < 0, looking for a possible fixed point. We have not found any stable solution, in 
agreement with the general arguments given in the introduction: the mixed point is indeed 
located in the region u > 0, and it is therefore irrelevant for the critical behavior of the 
dilute model. What remains to be checked is the stability of the O(M) fixed point. As we 
mentioned in the introduction, a general argument predicts that this fixed point is stable; 
the random (cubic) perturbation introduces only subleading corrections with exponent uj = 
—olm/vm- This exponent can be easily computed from 



10 



du 



(0,17), 



(3.27) 



which is ^-independent as expected. 



The analysis is identical to that performed for the stability of the Ising point in Ref. [92 



We use a conformal transformation and the large-order behavior of the series; then we fix 
the optimal values b opt and a opt of the parameters by requiring the estimates of v* and lj 
to be stable with respect to the order of the series used. The errors were obtained varying 
a and b in the intervals [a opt — 2, a opt + 2] and [6 op t — 3, b opt + 3]. As in Ref. [j9"2|| , the final 
result is reported with an uncertainty corresponding to two standard deviations. 

The final results for some values of M are reported in Table |Xil| , together with estimates 
of the theoretical prediction —olm/vm- F° r M > 3 these results clearly indicate that the 
0(M)-symmetric point is stable. The results are somewhat lower than the theoretical pre- 
diction, especially if we consider the high-temperature estimates of the critical exponents of 
Ref. [p. 12(1 . This is not surprising: indeed the estimates of the subleading exponents to show in 
many cases discrepancies with estimates obtained by using other methods. This is probably 
connected to the non-analyticity of the /9-function at the fixed point |M , 1 1 3 , 114 , |l 1 Q Jl 1 5|| . A 
similar discrepancy, although still well within a combined error bar, is observed for M — 2. 
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In this case, we obtain uj > 0, indicating that the fixed point is stable. The error however 
does not allow to exclude the opposite case. 
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